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"^ This is a continuation of Rcf.[l](arXiv:nlin. SI/0603008). In the present paper we review 

^H solutions to the modified Korteweg-de Vries equation in terms of Wronskians. The Wron- 

1-H skian entry vector needs to satisfy a matrix differential equation set which contains complex 

operation. This is different from the case of the Korteweg-de Vries equation. We introduce 

KH an auxiliary matrix to deal with the complex operation and then we are able to give complete 

C/j solution expressions for the matrix differential equation set. The obtained solutions to the 

(~| modified Korteweg-de Vries equation can simply be categorized by two types: solitons and 

•^ breathers, together with their limit cases. Besides, we give rational solutions to the mod- 

^ ified Korteweg-de Vries equation in Wromskian form. This is derived with the help of the 

Galilean transformed modified Korteweg-de Vries equation. Finally, typical dynamics of the 

p^J obtained solutions is analyzed and illustrated. We list out the obtained solutions and their 

t> corresponding basic Wronskian vectors in the conclusion part. 
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1 Introduction 

It is well known that the modified Kordeweg-de Vries (niKdV) equation, 

vt + 6ev'^Vx + Vxxx = 0, e = ±1, (1.1) 

played an important role in constructing infinitely many conservation laws [2] and Lax pair for 
the Korteweg-de Vries (KdV) equation. The Lax pair led to the breakthrough of the Inverse 
Scattering Transform (1ST) [3] and then the Soliton Theory. The mKdV equation is also fa- 
mous for its special soliton behavior, breathers. Actually, the mKdV equation, or its Galilean 
transformed version, 

Vt + 12evoVVx + 6eV^Vx + Vxxx = 0, e = ±1, (1.2) 

which is usually referred to as the mixed KdV-mKdV equation, arose in many physics contexts, 
such as anharmonic lattices[4], Alfven waves[5], ion acoustic solitons[6, 7, 8], traffic jam[9, 10], 
Schottky barrier transmission lines[ll], thin ocean jets[12, 13], internal waves[14, 15], heat pulses 
in solids[16], and so on. 

With regard to exact solutions, many classical solving methods, such as Hirota's bilinear 
method[17], the IST[18, 19], commutation methods[20] and Wronskian technique[21, 22, 23] has 
been used to solve the mKdV equation. For more references on solutions of the mKdV equation, 
one can refer to Ref. [20] and the references therein. 

In general, for a soliton equation with bilinear form, its solutions can be expressed through 
a Wronskian by imposing certain conditions on its Wronskian entry vector [24, 25, 26, 1]. For 
convenience in the following we refer to such conditions as condition equation set(CES). Usually 
for an (l+l)-dimensional soliton equation the crucial part in its CES is a coefficient matrix 
and the matrix and its any similar form leads to same solutions for the corresponding soliton 
equation. Thus it is possible to give a complete classification (or structure) for the solutions of 
the soliton equation by considering the canonical form of the coefficient matrix[26, 1]. It has 
been understood that the solutions generated from a Jordan form coefficient matrix are related 
to the solutions generated from a diagonal form coefficient matrix via some limiting proceure[l]. 
Therefore the latter solutions can be referred to as limit solutions. Actually, from the viewpoint 
of the 1ST, N solitons are identified by N distinct eigenvalues of the corresponding spectral 
problem, or in other words, A^ distinct simple poles {kj} of transparent coefficient ^jn^- When 
{kj} are multiple-poles, the related multiple-pole solution can be obtained through a limiting 
procedure like k2 — > fci from simple-pole solution. This limiting procedure is easily realized for 
solutions in Wronskian form[l]. Such a procedure is also helpful to understand the dynamics of 
limit solutions [27, 28, 29]. 

In Ref.[l] we mentioned four topics related to solutions in Wronskian form: to find the 
CES, to solve the CES, to describe relations between different kinds of solutions, and to discuss 
dynamics of the solutions. In the present paper, following this line, we will try to review the 
Wronskian solutions of the mKdV(e = 1) equation 

vt + Gv'^Vx + Vxxx = 0, (1.3) 

together with its Galilean transformed version, 

Vt + uvoVVx + ey Vx + Vxxx = o. (1.4) 

The main results of the paper are the following. 



• The CES for the mKdV equation can be given by 

Lpx =B(^, (1.5a) 

ift = - Aip^^j:, (1.5b) 

where (p is the A^-th order Wronskian entry vector, bar stands for complex conjugate and B 
is a nontrivial N x N constant complex matrix. We note that there is a complex conjugate 
involved in the CES and this makes difficulties when solving the CES. 

• There are only two kinds of solutions led by the CES (1.5): Solitons (together with their 
limit case) and breathers (together with their limit case). No rational solutions arise from 
(1.5) due to \M\ / 0. 

• Rational solutions to the mKdV equation (1.3) can be derived by means of the KdV-mKdV 
equation (1.4). 

• Dynamics of obtained solutions is analyzed and illustrated. 

The paper is organized as follows. In Sec. 2, we give a general CES of the niKdV equation 
and simplify the CES by introducing an auxiliary equation. Then in Sec. 3 we solve the CES 
and classify the solutions as solitons and breathers. In Sec. 4 we derive rational solutions to the 
mKdV equation. This is done with the help of the KdV-mKdV equation (1.4). Sec. 5 consists 
of dynamic analysis and illustrations. Finally, in the conclusion section we list out the obtained 
solutions and their corresponding basic Wronskian vectors. 

2 Wronskian solutions of the mKdV equation 

2.1 Preliminary 



An N X N Wronskian is defined as 



wi<Pu<p2,--- ,4>n) = \<PA^'\--- A^'^-'^l 



bf 4^) ... 4"-^^ 



^N 'Pn 'Pn 



(2.11 



where (/)■ = d'-cftj/dx and cj) = (0i, (f>2, ■ ■ ■ , (Pn)'^ is called the entry vector of the Wronskian. 
Usually we use the compact form [24] 

w{4>) = \4>, </.«,. . . , 4^-i)| = |0, 1, . . . , iV - 1| = |iV^l|, (2.2) 

where N — j indicates the set of consecutive columns 0, 1, . . . , A^ — j. In the paper we also 
employ the notation N — j to indicate the set of consecutive columns 1,2,... ,N — j. 

A Wronskian can provide simple forms for its derivatives and this advantage admits di- 
rect verification of solutions that are expressed in terms of Wronskians. The following matrix 
properties are usually necessary in Wronskian verification. 



Proposition 1. [30, 1] Suppose that 3 is an N x N matrix with column vector set {Sj}; f] is 
an N X N operator matrix with column vector set {^j} and each entry 0,j^s being an operator. 
Then we have 

N N 

where for any N -order column vectors Aj and Bj we define 



(2.3) 



A, o Bj 



{AijBij, A2JB2J, ■ ■ ■ , An,jBn,j) 



(2.4) 



and 






,^j_i, Ajor^j, ^j+i, 



,='N\- 



(2.5) 



Proposition 2. [24] Suppose that D is an N x {N — 2) matrix and a, b, c, d are Nth-order 
column vectors, then 



\D, a, b\\D, c, d\-\D, a, c\\D, b, d\-\-\D, a, d\\D, b, c\ 



0. 



(2.6) 



Besides, to write solutions of CES in simple forms one may make use of lower triangu- 
lar Toeplitz matrices. An A'^th-order lower triangular Toeplitz matrix means a matrix in the 
following form 
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NxN 



All such matrices form a commutative semigroup G]\f(C) with identity with respect to matrix 
multiplication and inverse, and the set Gn{C) = {A\ A G GAr(C), |^| / O} makes an Abelian 
group. Besides (2.7), we will also need the following block lower triangular Toeplitz matrix, 
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where A^ 



-Lj — , ^ J and {ttjs} are arbitrary complex numbers. All such block matrices 

also compose a commutative semigroup with identity which we denote by G^(C), and the set 
G^^{C) = {A^\ A^ G Gff{C), \A^\ / 0} makes an Abelian group, too. If all the elements 
are real, then we correspondingly denote the above mentioned matrix sets by G]\f(M.), GnO^), 
). For more properties of such matrices please refer to Ref. [1]. 



Gf^(M) and Gf^( 



2.2 CES of the mKdV equation 

By the transformation 

■ / , J \ ■ JxJ ~ J Jx /„ ri\ 

v = ^[ln-)^ = ^^J^, (2.9) 

the niKdV equation (1.3) can be bihnearized as[17] 

{Dt + Dl)f ■/ = {), (2.10a) 

Dlf-f = 0, (2.10b) 

where i is the imaginary unit, / is the complex conjugate of /, and D is the weh-known Hirota's 
bihnear operator defined by[31, 32] 

am an 

DTDla{t,x) ■b{t,x) = -—:^—a{t + s,x + y)b{t-s,x-y)\s=o,y=o, m,n = 0,l,--- . 

The bihnear mKdV equation (2.10) admits a solution / in Wronskian form. 
Theorem 1. A Wronskian solution to the bilinear niKdV equation (2.10) is given as 

f = W{<P) = \N^l\, (2.11) 

provided that its entry vector (p satisfies 

(t>x =B{t)4>, (2.12a) 

cl)t = -^xxx + C{t)ct>, (2.12b) 

where B{t) = {Bij(t))NxN andC{t) = {Cij{t))NxN OLfc twoNxN matrices oft but independent 
of X, and satisfy 

\B{t)\t = 0, tr(C(t)) GR(t), (2.13a) 

Bt{t) + B{t)C{t) = C{t)B{t). (2.13b) 

The proof is given in A. 

2.3 Simphfication of the CES (2.12) 

To solve the CES (2.12) with arbitrary B{t) and C{t) which satisfy (2.13) we first introduce a 
non-singular N x N complex matrix H{t) G C]\fxNit) such that ([33], also see [1]) 

Ht{t) = -H{t)C{t). (2.14) 

By H{t) we then introduce a new Wronskian entry vector 

^ = H{t)(k, (2.15) 

which transfers the CES (2.12) to the following, 

i^x =Bi>, (2.16a) 

V'i = -4Vxxz, (2.16b) 
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ifxx =A(/3, 




'Px =IB(^, 




(ft = - 4:ipxxX, 


'-^BP and 





where B = H{t)B{t)H~^{t) has to be a constant matrix independent of both x and t due to 
the compatibihty condition (2.13b) (noting that now C{t) = 0). 

We note that the Wronskians composed by (j) and tp, which we respectively denote by f{(j)) 
and f{ip), yield same solutions to the niKdV equation through the transformation (2.9) due to 
f{ip) = \H{t)\f{(j)). That means in the following one only needs to focus on the CES (2.16). 
However, since there exists a complex operation in (2.16), solutions can not be classified in 
terms of the canonical form of B, as done in [26, 1]. To overcome the difficulty we introduce an 
auxiliary equation 

V'xx = IV', (2.17) 

where A = BB. Our plan now is to first solve the equation set composed by (2.17) and (2.16b), 
which is nothing but the CES of the KdV equation and have been well studied in Refs.[26] and 
[1]. Then in the second step we impose condition (2.16a) on the obtained solution ip and finally 
get solutions for (2.16). 

In addition, noting that A and its any similar form A = P^^AP always generates same 
solutions to the mKdV equation, in the following let us focus on the CES 

(2.18a) 
(2.18b) 
(2.18c) 

where (/? = P~^'il), 

A = BB. (2.19) 

In the next section we will see that solutions to the mKdV equation can be classified in terms 
of the canonical form of A (rather than the canonical form of B). 

3 Solutions of the mKdV equation 

In this section, we list several possible choices of A and derive the related Wronskian entry 
vectors, one of which is for breathers. We will also discuss the limiting relationship of some 
solutions. 

First, for the matrix A defined by (2.19), we have the following result on its eigenvalues. 

Proposition 3. The eigenvalues of A defined by (2.19) are either real or, if there are some 
complex ones, appear as conjugate pairs. 

We leave the proof in B. 

According to the eigenvalues of A, we can categorize solutions to the mKdV equation as soli- 
tons and breathers, which correspond to real eigenvalues and complex eigenvalues of conjugate 
pairs, respectively. 

3.1 Solitons 

Case I. Solitons: When A has N distinct real positive eigenvalues {A^} its canonical form 
reads 

A = Diag(A?,Ai, •••, A^), (3.1) 



where, for convenient to discuss, we let Xj = £j\\kj\\ ^ 0, in which Ej = ±1 and kj can either 
real or complex numbers with distinct absolute values. We consider two subcases. 
(1). kj G M, i.e. 

A = Bmg{klkl ••• , k%). (3.2) 

Following the relation (2.19), we can take 

B = Diag(±A;i,±A;2, •••, ±kN). (3.3) 

We neglect the sign it because this can be compensated by the arbitrariness of kj. So next we 
take 

IB = Biag{ki,k2, • • • , /cat), kj € M. (3.4) 

For the above matrix B, the solution to the CES (2.18) can be 

ip= {Lpi,ip2,--- ,'fN)'^, (3.5a) 

where 

^j = aje^^ +iaje~^\ ^j = kjX - Ak^t + ^f , aj,aj ,kj,^f €R. (3.5b) 

When a+ = (-l)-''"\ aj = 1, (3.5b) reads 

cp. = (-iy''^e^^+ie-^^. (3.6) 

In this case, the corresponding Wronskian solution (2.9) with /(</?) can be written as[34] 

where the sum over ^u = 0, 1 refers to each of fij = 0, 1 for j = 1, 2, • • • , iV, and 

This coincides with the A'^-soliton solution in Hirota's exponential polynomial form[37]. 
(2). kj = kji + ikj2 G C, i.e. 

A = Diag(A?, Xl • • • , A^), X] = fcfi + k%. (3.7) 

In this case we have 

B = Diag(A;i, A;2, • • • , /cat), kj e C, (3.8) 

and the solution to the CES (2.18) can be given by (3.5a) with 

^j = ij (a+e«^- + iaj e'^^ ) , ^j = XjX - AX]t + if , a+ , a" , (,f G M, (3.9a) 

where 

7, = 1 + ^^^V^' A,=e,P,||. (3.9b) 



k 



i2 



Now, if we compare (3.5b) and (3.9a), we can find that both of them provide same solution to 
the mKdV equation through the transformation (2.9). In particular, when kji = 0, i.e., 



B = Diag(iA;i2,iA:22, ■■■,ikN2), 
and we take Xj = kj2, the entry function (3.9a) reduces to 



(fj 



[1 + i)(a+e«^- + ia-e-^^), ^^ = kjX - Akh + a'\ a+ oT ,£!f^ G R. 



,(0) 



J '^J 



(3.10) 



(3.11) 



This was first given by Nimmo and Freeman[22] as Wronskian entries for soliton solutions. 
Let us remark this case as follows. 

Remarks: 

• WhenM in the CES (2.18b) and (2.18c) is a diagonal matrix (3.8) of which the diagonal 
elements have different absolute values, no matter these diagonal elements are real or 
complex, the related Wronskian f{ip) generates N -soliton solutions to the niKdV equation. 



Besides, we specify the non- degenerate condition* 



(3.12) 



i.e., ki and kj (i ^ j) can not appear on the same circle with the original point as the 
center of the circle. Otherwise, the solution degenerates. 



Case II. Limit solutions of solitons: Corresponding to Case I, we discuss two subcases. 
(1). Let 

/ kf O--- 0\ 



A 



2ki kf 
1 2ki kf 








, fci G M. 



(3.13) 



1 2ki kl J 



N 



V 

In this subcase, general solution to the equation set (2.18a) and (2.18c) is[l] 

where 

1 






(3.14) 



(3.15) 



^1 is defined in (3.5b), Gn{C) is the commutative set of all the A^th-order lower triangular 
Toeplitz matrices, see Sec. 2.1. 



*The above non-degenerate relation can also be described as follows. Define the equivalent relation ~ on the 
complex plane C: 



kj iff 1 1 fci 1 1 = j I fej I 



The quotient space C/ ~ denotes the positive half real axis. Then to get a non-degenerate soliton solution one 
needs ki '^ kj, {i ^ j). 



For the matrix B satisfying (2.19) we take 



B = 


/ fci 
1 ki 


••• \ 
••• 


V 


••• 1 ki ] 


Then, substituting into (2.18b) yields 




AQl^+BQ,^, 


= M{AQ+ + BQ^). 


Meanwhile, it can be verified that 






Q-t. = mt 


> Qo,x = -®2o"- 



(3.16) 



N 



(3.17) 



(3.1J 



Noting that B G Gn(^)-, substituting (3.18) into (3.17), and making use of the commutative 
property of G]\f{C), we have 



B(^Q+ - BQq) = B(^Q+ + BQq 
Then, compared with (3.17) we immediately get 

A = A, B = -B, 



(3.19) 



(3.20) 



which means A is real and B pure imaginary. In the end, the solution to the CES (2.18) can be 
described as 

if = A+Q^ + iA-Qo, ^^gGjv(M). (3.21) 

We note that solutions generated from (3.21) can also be derived from the solution given in 
Case 1.(1) by a limiting procedure (cf.[l]). Let us explain this procedure by starting from the 
following Wronskian 



fiv) 



W{ipi,(p2,--- ,'^n) 



iN 



Yr=2ikj - hy-' 



(3.22) 



with c/?i = ipi{ki,x,t) = afe^^ +ia]^e^^^ as defined in (3.5b) and ipj = (pi{kj,x,t). (3.22) gives a 
Wronskian solution to the bilinear mKdV equation (2.10). Taking the limit kj — )■ ki successively 
for J = 2, 3, • • • ,N and using L'Hospital rule, the Wronskian (3.22) goes to a Wronskian W{f) 
with entry vector (3.21), where the arbitrary coefficient matrices A,B can come from a^ by 
considering a^ to be some polynomials of ki (cf.[l]). 

(2). Corresponding to Case 1.(2), let us consider 



/ kl,+k^ 



A 



12 

2ku ^=^ 





'^11 ^ '^'12 
12 








1 2^2 kl^ + kl^ 
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1 2ki2 kf^ + kJ2 J 































(3.23) 



N 
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where ku, ki2 G M and ki2 7^ 0. It then follows from (2.19) that one can take 



/ -iki 
1 -iki 



V 







\ 



1 -iki J 



, ki = fell +iki2. 



(3.24) 



N 



For the matrix A defined by (3.23), the general solution to the equation set (2.18a) and 
(2.18c) can be written as 

^ = ^p+ + 6p-, A^BeGNiC), (3.25) 

where^ 



V'. 



1 



V 0,0' '0,1' ^ ' 0,N-l) ' ' 0,s — g]'^ki2'^ ' 



6 =Xix- 4Xit + er , Ai = V ^11 + k 



'-12- 



(3.26a) 
(3.26b) 



Next, we substitute (3.25) together with (3.26) into (2.18b) so that we identify A,13 for (2.18b). 
This substitution yields 



AP+, + BV^^, = M{Ar+ + BV^) = ABV^ + BMV^ , 
where we have made use of the commutative property of G]\f{C). Then, noting that 

where 





v+ 


= WV+, 


^0,x 
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-WVo 
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(3.27) 



(3.28) 



W 

\ WN-l WN-^2 WN-3 ■■■ Wl Wo J j^ 

and substituting (3.28) into (3.27), we get the relation 

Aw = M, -Bw = m. 

To have a clearer result, we write 

^ = ^1 + iA2, B = Bi + iB2, IB = Bi + ilB2, 
where Ai,A2,Bi,B2 are in GnO^) and 

\ 



1 



^...^1' 



(3.29) 



(3.30) 



(3.31) 



/ ki2 

1 ki2 

\ 



1 ki2 ) 



( -ku 

-ku 

V 



\ 



(3.32) 



-ku J 



^Here we define Vq^ by taking derivative with respect to fci2 rather than fcn because in this case we always 
have fci2 7^ 0. 
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Then (3.30) yields 

Ai{W - Bi) - A2M2 = 0, A1B2 - ^2(Bi + W) = 0, (3.33) 

l3i{W + Bi) + S2B2 = 0, B1M2 - ^2(Bi -W) = 0, (3.34) 

which provides the equation sets for Ai,A2, Si, ^2- We note that since all the elements in (3.33) 
and (3.34) are in G]\f(R), we can treat (3.33) and (3.34) as ordinary linear equation sets. Let us 
look at (3.33). It first indicates 

A2=M^^{W-Mi)Ai, (3.35) 

and further, by eliminating A2, 

We note that 

In fact, by calculation we find 



ll + Ml-W'^)Ai = 0. (3.36) 

mI + mI-W"^ = 0. (3.37) 



^"^^''-{Tl^drjikl + kl,), ^>J, 

from which it is easy to verify (3.37). Thus, it turns out that ^1 can be an arbitrary element in 
GAr(M) and then 

A = AiM, M = In + M2^{W - Bi), (3.38) 

where In is the A^th-order unit matrix. Similarly, we can find 

B = iB2M, S2GGjv(M). 

Thus we can conclude that the solution to the CES (2.12) can be given by 

^ = 7W(^+pQ+ + i A-Vq), A"^ G Giv(M), (3.39) 

where M is defined in (3.38). Obviously, the matrix M contributes nothing through the trans- 
formation (2.9) to the mKdV equation. Therefore in practice we may remove M and use the 
effective part 

^ = ^+p+ + i A-Vq, A"^ G G7v(M), (3.40) 

which is as same as (3.21). 

3.2 Breathers 

Case III. Breathers: When A has distinct complex (conjugate-pair) eigenvalues, we may have 
breather solutions. Let us consider a 2A'^-th order matrix, 

A = Diag(A;2, kl . . . , k%, k%)2N, kj / 0, (3.41) 

where kj £ C, j = 1,2,- ■ ■ ,N. The matrix B which generates breathers is 

B = Diag(ei, 62, • • • , eN)2N, e, = ( ^ ^^' ] . (3.42) 

12 



For this case, the Wronskian entry vector, i.e., the solution to the CES (2.18), can be taken as 

^ = (V'll) 9512, V'21, V'22, • • • , '^Nl, fN2)'^, (3.43a) 

where 






■ A 



^j^ + bjC ^\ (pj2 = CLje^^ — bjS 



..?. 



,(0) 



kjX-Akn + ^Y', aj,bj,a' ^C. 



(0) 



"3" ' ^J 



(3.43b) 



With the above 93 as basic entry vector, the Wronskian f{(p) wih provide breather solutions for 
the niKdV equation. For non-trivial solutions we need kji 7^ 0, (see Sec. 5. 3). 

Case IV. Limit solutions of breathers: In this case let us consider the following block 
matrix, 

//C 0... o\ 



A 



h K, K, 



\ 








where 



'==1 S° !• ^ 



2ki 
2ki 



h K, IC )^^ 
h-- 



1 
1 



The matrix B satisfying (2.19) can be taken as 



I B ... \ 
I2 B ... 

\ ... I2 B J 



2N 



where 



1 



, B 



1 

The Wronskian vector of this case is also in the form of 



fci 
Ail 



V' = (V'l,!) 9^1,2, 9^2,1, 9^2,2, • • • , ^N,l,^N,2y 



and for convenience we set 

<f^ = (V'l,!, 9^2,1, • • • , VN,i)'^, f~ = (9^1,2, 9^2,: 



,^N,2) 



Substituting (3.46) into (2.18) with the above A and B we get 



^t 



^Nf , fa: 



mf 



ft 



^fxxx^ 



(3.44a) 



(3.44b) 



(3.45a) 



(3.45b) 
(3.46a) 

(3.46b) 

(3.47a) 

(3.47b) 
(3.47c) 
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where 



An 



f kf 

2ki kl 

1 2ki kl 



V 



/ kl 
1 





kl 



V 













\ 







(3.48a) 



1 2ki k-f I ^ 

\ 




(3.48b) 



1 kl ) 



N 



From (3.47a), (3.47c) and (3.48a), one first has 

cp+ = AQ++BQo, ^~ =CQ++VQ^, A,B,C,V £ GNiC), (3.49) 

where Qq are defined by (3.15) together with ,^i defined in (3.43b). Then substituting (3.49) 
into (3.47b) and making using of (3.18), one finds 

C = A, V = -B. (3.50) 

Thus, (3.49) reads 

^+ = AQt + BQ^ , i^- = AQ+ - BQ^ , A^gGjv(C). (3.51) 

The Wronskian f{^) with 99 defined by (3.46) with (3.51) will provide a limit solution of 
breathers for the mKdV equation. In fact, quite similar to the procedure for the limit solutions of 
solitons we described in Sec. 3.1, here the Wronskian /(</?) with (3.46) is related to the A^-breather 
solution by taking the limit kj — )■ fci successively for j = 2, 3, • • • ^N . 

4 Rational solutions 

4.1 Backgrounds 

Following solution structures of the KdV equation[l], rational solutions should be led from zero 
eigenvalues of the coefficient matrix A. For the mKdV equation this requires |B| = 0. However, 
this is not allowed in the Wronskian verification (see Appendix C.4) because a trivial B will 
lead to a zero Wronskian in the light of (A.l). Thus, it is clear that to get non-trivial rational 
solutions to the mKdV equation, we need a non-trivial matrix B. 
Let us go back to the KdV-mKdV equation (1.4), i.e., 

Vt + l2vQVVx + QV^Vx + Vxxx = ^, (4.1) 

which is related to the mKdV equation (1.3) through a Galilean transformation 

v{x,t)=VQ + V{X,t), x = X + QvQ^t, (4.2) 

where vq is a real parameter. We note that the equation (4.1) admits non-trivial and non- 
singular rational solutions when vq 7^ 0. Then using the transformation (4.2) rational so- 
lutions to the mKdV equation can be obtained. This fact has been realized via Backlund 
transformation(BT) [35] and Hirota method with a limiting procedure[36, 37], but the presen- 
tation for high order rational solutions is complicated. In the following we will derive rational 
solutions in terms of Wronskian, which provides not only explicit but also impact forms for high 
order rational solutions. 
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4.2 Rational solutions 

Still employing the same transformation as (2.9), i.e., 

y = i(ln^)^, (4.3) 

the KdV-niKdV equation (4.1) can be written into the bilinear form[38, 39] 

(A + I?i)/-/ = 0, (4.4a) 

{Dl-2ivoDx)f-f = 0. (4.4b) 

For the solutions to (4.4) in Wronskian form, we have 
Theorem 2. The bilinear equation (4.4) admits Wronskian solution 

f = W{cl)) = \N^l\, (4.5) 

where the entry vector satisfies 

i(px = vo(j) + B{t)4>, (4.6a) 

(pt = -40XXX + C{t)ct>, (4.6b) 

in which B{t) and C{t) are two N x N matrices oft but independent of x, and satisfy 

\Bit)\^0, (4.7a) 

tr(C(t)) G M(t), (4.7b) 

Bt{t) + B{t)C{t) - C{t)B{t) = 0. (4.7c) 

The proof is similar to the one for Theorem 1, but (4.6a) results in complicated expression 
for / and its derivatives. We leave the proof in D. 

As in Sec. 2. 3, with the help of the auxiliary matrix A and auxiliary equation 



ipxx = A99, (4.8) 



one can simplify the CES (4.6) to 



iifx = VQip + By?, (4.9a) 

(ft = -4:ipxxx, (4.9b) 

where both A and B are N x N complex constant matrices and related by 

A = BB-t;^/Ar, (4.10) 

in which Ij^ is the A'^th-order unit matrix. It might be possible that here we discuss all possible 
solutions according to the eigenvalues of A, as we have done in the previous section. However, 
since we have had a clear description for solitons and breathers of the mKdV equation in the 
previous section, and the parameter vq will bring more complexity, in the following we can 
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neglect the discussion of A and let us only focus on rational solutions. In fact, rational solutions 
correspond to the zero eigenvalues of A. 

We derive rational solutions as the limit solutions of solitons. The A^-soliton solution corre- 
sponds to 

B = diag( - ^vl + kl -yjvl + kl--- ,-^vl + kl), (4.11) 

where we take kj to be A^ distinct real positive numbers. In this case, a solution to the CES 
(4.9) is 

V = (v3i,V2,--- ,VnY , (4.12) 

with 

Lpj = ^j2vQ + 2ikj e'l' + y^2wo - 2ikj e"''^ , r]j = kjX - Ak^t, (4.13) 

This provides an A^-soliton solution to the KdV-mKdV equation (4.1) through the transforma- 
tion (4.3) with /((/>). 

Now in the CES (4.9a) we take B to be a lower triangular Toeplitz matrix, 

/ ao ^ 



with 



a, 



\ aN-i ■■■ a\ OH) j 

Jvl + kl , , (i = 0,l,---,iV-l). 

V k■^ ^0 



fei=0 



1 d^^ 

''- {2j)ldki'^ 
In this case the CES (4.9) admits a solution 

V' = (V'i,V'2,--- ,i^NV 

with 



^ 



'i+i 



J2jy.dk^ 



"pi 



ij = 0, 1, 



,N-1) 



(4.14) 

(4.15) 

(4.16a) 
(4.16b) 



fci=0 



where (pi is defined by (4.13). Then, with such a ^ as basic column vector, the Wronskian 

/ = VF(V') = \N^l\ (4.17) 

provides non-singular rational solutions to the KdV-mKdV equation (4.1). A simplified form of 
these solutions is 

-2(Fi,xi^2 - FiF2,x) 



V{X,t) 



F? + Fi 



Fi=Re[f], F2=lm[f]. 



We list the first three non-trivial / (for A^ = 2, 3, 4, respectively) in the following, 
f=4{2voX + i), 



f 
f 



16i;o-v/2t;o 
3 
1 4X2 



X-^ + 12t 
16 



4«2 + 2t;o^ ^4vl 



Vq 4 3 



{l2tX + X*) - —v^ (720*2 - 60tX^ - X^) 



45 



+ 4ron 



24t X 



16 



+ -r+32tX' + —Xn. 



% % 



15 



(4.18) 

(4.19a) 
(4.19b) 

(4.19c) 



For the rational solutions to the mKdV equation (1.3), we have 
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Proposition 4. The non-trivial rational solutions to the niKdV equation is given by 

v{x, t) = vo ^2 I j72 ' X = x- ^Vot, (4.20) 

where f is the Wronskian (4.17) composed with (4.16), Fi = Re[/], F2 = Ini[/]. 
The first non-trivial rational solution to the niKdV equation reads 

^ = ^o-^^P^, X = x-6vlt, (4.21) 

and the second one 



vo- , ,,.., , .J%v,, !'' ., , 1 ,, , X = x-Qvlt. (4.22) 



2v^ l&vf^ 

4^2(^3 + 12t - 1^)2 + 9(X2 + 1)2 ' 



We note that there is a bilinear BT[22] related to (4.4) and the BT admits multi-soliton 
solutions[22] and rational solutions[40] in Wronskian form. Besides, rational solutions of the 
mKdV equation (for both e = ±1) may also be derived (through the Galilean transformed 
version) by using the long-wave-limit approach described in Ref.[41] (also see [23, 42]). 

In next section we will discuss dynamics of solutions including these rational solutions. 

5 Dynamics analysis 

In this section we investigate dynamics of two-soliton solutions, limit solutions, breathers and 
rational solutions. To describe the relationship between two-soliton solution and the simplest 
limit solution, let us start from the asymptotic behaviors of two-soliton interactions. 

5.1 Dynamics of solitons 

Soliton solutions to the niKdV equation (1.3) can be described by 

t'(x,t) = -^^^m^|P^, Fi=Re[/], F2=Im[/], (5.1) 

where / = f{ip) = |A^ — 1| is the Wronskian composed by the basic column vector (p which is 
defined by (3.5), or equivalently, by either (3.9a) or (3.11) in Case I. We note that from the 
transformation (2.9) the solution to the mKdV equation (1.3) can also be written as 

t; = -2(arctan-i) = 2(arctan -^) , Fi = Re[/], F2 = Im[/], (5.2) 

V F2 ' X V F\ J X 

while (5.1) gives a more explicit form. 

In the following we investigate 1- and 2-soliton solutions, which are, respectively, correspond- 
ing to 

f = fi=f{^i) = ^1, (5.3a) 

f = h=m^i,^2f)= I' l'^^ , (5.3b) 

9^2 ^2,x 
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with ifj defined by (3.5b), i.e., 



^j = a+e«^ + ia^ e'^^ , ^^ = %x - Ak^^t + q" , a^ , k^ , Cf ' G 



Then, one-soliton solution to the mKdV equation (1.3) reads 



,(0) 



-2 • sgn 



• ki ■ sech( 2A;ix - Sk\t - In 



(5.4) 



(5.5) 



as depicted in Fig.l, where for convenience, we call (a) soliton and (b) anti-soliton due to the 
signs of their amplitudes. 





(a) (b) 

Figure 1: The shape of one-soliton solution given by (5.5) at t = 0. (a). Soliton for oJi'' = a]~ = 1, 



ki = —1 and C| = 0. (b). Anti-soliton for a^ 



1, ki = 1 and ^\ 



(0) 



0. 



Obviously, this soliton is identified by the amplitude 

Amp = —2 • sgnT^l • ki 

and top trace (trajectory) 

x{t)=Aklt+--\n\'^\, 
2ki ' a|f ' 

or velocity Akf. Obviously, solitons of the mKdV equation are single-direction waves. 

Next, let us look at two-soliton solution. The two-soliton solution of the mKdV (1.3) can be 
expressed by (5.1) where from (5.3b) 



(5.6) 



Fi = {k2 - A:i)(ara2-e4('=?+'=')*-(^i+'=2)" + a+a+e-4(*^?+*^')*+(*^i+'=2)"), 



(5.7a) 



F2 = ik2 + /fci) (ara+e^('=?-'=2)*-('=i-'=2)" - a+a^e^(-'=?+'=2)*+('=i-'=2)-) . (5.7b) 

We assume that a- 7^ and set sgn[^] = sgn[^], so that Fi,F2 are not zero at same time. 

Without loss of generality, we also take < \k2\ < \ki\. For the analysis of asymptotic behaviors, 
it is convenient to use the following expression 



V = 2{ arctan 



F2 
FiJx 



(5.8) 
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(a) (b) 

Figure 2: Two-soliton interactions, (a). Soliton-soliton for a^ = 



a2 



(0) 



ki = 0.8, k2 = 0.5 and Ci 
ki = -0.8, /C2 = 0.5 and Ci 



(0) 



S2 



'?2 



0. (b). Soliton-anti-soliton for a\^ 
= 0. 



02 



1, &2 = 

= 1, 62 



There are two types of 2-soliton interactions, soliton-soliton (or anti-soliton-anti-soliton) 
interaction and soliton-anti-soliton interaction, as shown in Fig. 2. 

To investigate asymptotic behaviors of the two solitons involved in interaction, we first name 
them fci-soliton and A:2-soliton, respectively. Then we rewrite the two-soliton solution (5.8) in 
terms of the following coordinates, 



[Xi = x-4klt,t), 



which then gives 



a^ a2 e 



2kiXi 



a> 



+ ^-g8fc2(fef-fci)t+2fc2Xi)(ju^ + fca) 






(5.9) 



(5.10) 



Noting that for any < 1^2 1 < \ki\ it is always valid that /cj — fc^ > and ^r4fc^ > 0, we can 
keep Xi to be constant and let t go to infinity. Then we can find there is only /ci-soliton left 
along the line Xi = const, and also find how the fci-soliton is asymptotically identified by its 
top trace and amplitude, for both t — )• ±00. 

As for details, when ^2 > 0, t — )• —00 or A;2 < 0, t — )• +00, i.e. sgn[/i;2] • t 
(5.10) becomes 



-00, the solution 



2kiXi 



V =2( arctan 



-2 ■ sgn 



af{ki + /c2)e 

a-iiki — k2) ^ ^1 

ki • sech( 2A;iXi — In 



+ ln 



ki + /C2 
ki - /C2 



(5.11) 



and when ^2 > 0, t — )■ +00 or A;2 < 0, t — >• —00, i.e. sgn[A:2] • t — ;■ +00, (5.10) becomes 



-a^ {ki + ^2) 



V =2[ arctan , „, ,^, , , 

V ate^^^^{ki-k2)^x^ 



-2 ■ sgn 



fli 



ki ■ sech( 2kiXi — In 



In 



fcl + /C2 

ki — k2 



(5.12) 
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We can also rewrite the two-soliton solution (5.8) in terms of the coordinates 

{X2 = x-Aklt,t), 
and do a similar asymptotic analysis for the /c2-soliton. Finally, we reach to 



(5.13) 



Theorem 3. Suppose that sgnT^l = sgnT^l, a- / and < \k2\ < \ki\ in (5.10). Then, 
when sgn[A;2] • t — t- ±00, the ki-soliton asymptotically follows 



top trace : x{t) = Ak^t H — — In 



, 1 , ki+k2 

± —— m — , 

2ki Ki — k2 



amplitude : 2 • sgn 



• ki, 



(5.14a) 
(5.14b) 



and when sgn[/ci] • t — )■ ±00, the k2-soliton asymptotically follows 



top trace : x(t) = 4^2^ H r~ ^^ 

2k2 



amplitude : — 2 • sgn 



fco. 



^ 1 , ki + k2 
± — — m — . 

2fc2 Kl - K2 ' 



1 Iv, k2+ki 



(5.15a) 
(5.15b) 



The phase shift for the kj-soliton after interactions is ^ ^^ k-k • 



Now it is completely clear how the two-soliton interactions are related to the parameters 
{fcj, Oj, 6j}. This will be helpful to understand the asymptotic behavior of limit solutions. 



5.2 Asymptotic behavior of limit solutions 

The simplest limit solution in Case II is 

a^a^; (-48A:ft + 4A;ix) 



V = 2{ arctan ■ 



^-2g8fcft-2fcia; , ^+'^ ^-8kft+2ki 



x/ X 



(5.16) 



This is derived from (5.2) with 



where ipi is defined in (3.5b), i.e.. 



-,+J 



^1 ^l,x 

dki^i {dki^i)x 



(0) 



,(0) 



ip^ = a+e''^ +ia^e~''\ ^^ = k^x - Akp + C\^' , af,ki,^\'''€ 



(5.17a) 



(5.17b) 



The solution (5.16) is depicted in Fig. 3. We characterize dynamics of the solution by the 
following two points, which are typically different from the interaction of two normal solitons 
that we described in the previous subsection. These two points are 

• Soliton-anti-soliton interaction with (asymptotically) asymmetric wave shape, 

• Top trace of each soliton is asymptotically governed by logarithm and linear functions. 
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ISvN. 155^ IfiJ l>6i 17^J 17: 



(b) 



Figure 3: Limit solution given by (5.16) for af = a^ 
motion, (b). Asymmetric wave shape at i = 160. 



,(0) 



1, ki = 0.5 and ^J = 0. (a). Shape and 



,± 



The first point can be explained as follows. Recall the two-soliton interaction with Og 



a-|^ 7^ and ki-k2 > 0. According to Theorem 3, this is a soliton-anti-soliton interaction and the 
absolute of amplitude of each soliton is 2\kj\. Obviously, the asymmetric shape of limit solution 
(5.16) coincides well with the limit k2 — t- fei. 

To understand the second point, again we put the limit solution in the coordinate system 



and this gives 

V = 2( arctan 

which is described in Fig. 4. 



Y = x-Aklt,t), 
Aafa^e^''^^ki{Y -8kft) 



(5.18) 
(5.19) 






10 1^ 



Figure 4: Limit solution given by (5.19) for a'^ 
(b). Asymmetric wave shape at t = 160. 



(b) 



1 and ki = 0.5. (a). Shape and motion. 



For convenience here we suppose that A;i > 0. Then by analyzing the leading terms as 
t — )• ibcxD in the numerator and denominator in (5.19), we can conclude the asymptotic behaviors 
of the limit solution (5.19) as follows. 
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Theorem 4. Suppose that a^ ^ and ki > in (5.19). When t — t- — oo there are two waves 
moving with amplitudes =f2 • sgn[^l • ki and top traces governed by the logarithm functions 



Kt = ^ In 
^ 2ki 



± ln{-32kft] 



(5.20a) 



where the subscript ^ of Y stands for Y — )■ ±00. When t — )■ +00, there are also two waves 
moving with amplitudes ±2 • sgn[— ^] • fci and top traces 



Y± 



2/ci 



In 



± ln(32/cf t) ) . 



(5.20b) 



Replacing y by x by using (5.18) the top traces of the waves in Fig. 3 are then asymptotically 
governed by linear and logarithm functions of t. 

Let us end up this subsection with the following remark. With regard to the limit solutions 
the top traces (or waves trajectories) are governed by (linear and) logarithm functions should 
be a typical characteristic, which differs from normal soliton interactions with straight line 
trajectories. (See [30, 27, 28, 29] for more examples). 

5.3 Breathers 

Wronskian entries in Case III provides breather solutions to the mKdV equation (1.3). The 
simplest one corresponds to 



where 



/ 



9^11 
9^12 



^12,2 



Fi + iF2 



V^ii = aie^^ + bie ", ipi2 = aie" — bie ", 
^1 = kix — Akit, ki,ai,bi € C, 



(5.21a) 



(5.21b) 
(5.21c) 



Fi =4A;ii(aii6ii + 012612) cos(24A;iifci2t - 8/012* - 2A;i2x) 
+ 4^11(012611 - 011612) sin(24/ciifci2t - 8kf2t - 2ki2x) 

F2=- 2^126-2^=11 (4'=?i*+12fc?2t+-) [(52^ + 62^)el6fcn* 
111 + ^12) 
and we have written 



+ (o2 +o2)e4'="(i2fc2^*+-)l 



(5.21d) 
(5.21e) 



h = kn + iki2, oi = Oil + iai2, h = 611 + i6i2. 
Then the breather solution is expressed by (5.1) with the above Fi, F2. We further assume that 
011611 + 012612 



sin( 



a 



a 



\/(aii6ii + 012612)2 + (012611 - 011612) 
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and then rewrite (5.1) as"'' 



-2( arctan — ) , 



where 



P = 2kna sin [2ki2{x - 4t{3kf^ - /c^s)) - ^] > 

Q = ki2[{ali + a22)e2'="(-+4*(3'=f2-fc?i)) + (ftf^ + i,2^y-2k,,ix+it(3q,-kl,))y 

Such a breather is described in Fig. 5. 



(5.22a) 

(5.22b) 

(5.22c) 



^-.>~^.i\l 





(c) 

Figure 5: Breathers given by (5.22). (a). A moving breather with ai = 1 + i, 6i = 2 + 0.5i, 
ki = 0.8 — 0.6i. (b). A stationary breather with oi = 6i = 1, ki = v3 + ^- (c). A 2D-plot of (b) 
at t = 0. 

Fig. 5 (a) shows an oscillating wave moving along a straight line. The oscillation comes from 
the sine function and the frequency depends on both x and t. To understand more on the wave 
we use the coordinates 

{Z = x + At{3kl2-kl^),t) (5.23) 



to rewrite the solution (5.22) as 



2 1 arctan 



2/ciiasin {2ki2{Z - 8t{kl^ + kf^)) 



(5.24) 



ku ((afi + af2)e2'="^ + [bf^ + bl2)e-^k^,z) ; ^' 
and then fix Z = 0, i.e., looking at the wave along the straight line Z = 0. Then it is clear that 

• The breather travels along the straight line Z = 0, in other words, the wave speed is 
4:{3kf2 — kfi), which means it admits bi-direction traveling. 

• A stationary breather appears when 3A;f2 = ^ii? as depicted in Fig. 5(b). 

• Under the coordinate system (5.23) the frequency only depends on t and the period reads 



T 



IT 



Skuikfi + kf^) 



(5.25) 



*Froin (5.22) we can see that kij 7^ is necessary for getting nontrivial breather solutions. 
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We note that in two dimensions (fixing t) the breather is in fact a spindle-shape wave. Let 
us go back to the solution (5.22). If we fix time t, then the breather oscillates with frequency 
L^ and its amplitude decays by the rate e~'^'^'^^^' as x — t- ±00. That means if |A;i2| is small and 
I fell I is relatively large so that the amplitude decay becomes the dominating factor, we get a 
'normal' breather as shown in Fig. 6(a); while if |A;ii| is small enough and |A:i2| >> |A;ii| so that 
the oscillation dominates, we will see a spindle-like wave shown in Fig. 6(b). In the latter case, 
the wave will travel with high speed and high oscillating frequency. This can lead to overlaps of 
'normal' oscillating waves (like Fig.6(a)) during their traveling, which makes a spindle shape. 





(b) 

Figure 6: Shape of breathers given by (5.22). (a). 'Normal' shape at t = with ai = 1 + i, 61 = 
2 + 0.5i, ki = 0.8 — O.Qi. (b). Spindle-shape breather at t = with oi = 1 + i, bi = I + i, 
ki = 0.5 + 5i. 

Finally, in this subsection we list and depict two-breather solution and the simplest limit 
breather solution, without further asymptotic analysis. Both solutions can be given by (5.1) 
with / being a 4 by 4 Wronskian 






where for the 2-breather solution 



.ii 



ifji = ttjC^^ + bjC ^\ (pj2 = CLje^^ — bjC ^\ ^j = kjX — Akjt + C 
in which Cj, bj, kj,^- S C, and for the limit breather 

if = (ifii, ipi2, dk^Lpii,di^Lpi2)'^ , 



(0) 



(5.26) 

(5.27a) 
(5.27b) 

(5.28) 



in which ipii^ipi2 are defined by (5.27b). 

Fig. 7 shows the two-breather interaction where from the density plot (b) one can clearly 
see that the two breathers are traveling along straight lines and a phase shift appears after 
interaction. Fig. 8 shows the shape and motion of a limit breather solution, where from the 
density plot (b) one can clearly see that the breather trajectories are not any longer straight 
lines. Here we conjecture that they are governed by logarithm functions. 
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(a) (b) 

Figure 7: Shape and motion of two breather solution given by (5.1) with (5.26) and (5.27). (a). 
3D-plot for ai = bi = a2 = b2 = 1, h = 1 + 0.5i, k2 = 0.8 - 0.6i and cf^ = ^f ^ = 0. (b). 
Density plot of (a) for x G [-12, 12], t e [-4,4]. 





Figure 8: Shape and motion of the limit breather solution given by (5.1) with (5.26) and (5.28). 



(a). 3D-plot for oi = 6i = 1, ki = 0.8 + 0.5i and ^[ 

xe [-10,10], te [-8,8]. 



(0) 



0. (b). Density plot of (a) for 



5.4 Dynamics of rational solutions 

The first non-trivial rational solution to the mKdV equation is (4.21), i.e., 

4fo 



V = vq 



Av^X^ + 1 ' 



X = X — QVnt. 



(5.29) 



This is a non-singular traveling wave moving with the constant speed 6wg, constant amplitude 
-3^0 and asymptotic line v = vq. It is depicted in Fig. 9. 
The next rational solution is given by (4.22), i.e.. 



V = vo 



UvoiX^ + AX^ 



3 



- 24Xt) 



4T;2(X3 + 12t-|4)2 + 9(X2 + 






, X = X — 6vQt. 



(5.30) 
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Figure 9: Shape and motion of the rational solution given by (5.29) for vq 
in (b) and uq = 1 in (c). 



It can be viewed as a double-traveling wave solution 



V = vo 



4t;2(x3-|ll)2 + 9(x2 + 






(5.31) 



with 



X = x- Gv^t, Y = x- 22vlt. 



(5.32) 



However, it does not show interactions of two single rational solutions. Only one wave is left for 
large x,t (see Fig. 10). We re-depict Fig. 10(a) in Fig. 11 by a density plot so that we can see the 
wave top trace clearer. 





(a) 



(b) 



Figure 10: Shape and motion of the rational solution given by (5.30) for vq = —0.8 in (a) and 
Vq = 0.8 in (b). 



To realize the asymptotic behavior analytically, we rewrite the solution (5.30) in the following 
coordinates system 



(X, T = X3 + 12t + — J-) 



IQv^X' 



(5.33) 
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(a) (b) 

Figure 11: (a): Density plot of Fig.lO (a) with X G [-10, 10], t G [-30,30]. (b): (a) overlapped 
by the trajectory curve given by (5.36). 



and this gives 



l2v,{2,X^ + i^X^-2TX + ^) 



Then, by calculation it can be found that for given large X the wave (5.34) has a single stationary 
point at T = where v gets a local extreme value 

_ vo{l - 12^2x2^ 



1 + 4i;2x2 



(5.35) 



which goes to —'ivQ as X — t- ±00. Thus we can conclude that for large X,t the wave asymptot- 
ically travels along the curve 

r = x3 + 12t + -l- = (5.36) 

IgVqX 

with amplitude — 3fo. Fig.ll(b) displays a density plot overlapped by the above wave trajectory 
curve. 

More details on the rational solutions to the niKdV equation can be found in [40] where the 
rational solutions are derived via bilinear Backlund transformation. 

6 Conclusions 

6.1 Summation 

In the paper we reviewed the Wronskian solutions to the mKdV equation (1.3) in terms of 
Wronskians. When a solution is expressed through the Wronskian 

/ = /(v5) = |iV^l|, (6.1) 

one needs to solve the finalized CES (2.18) together with (2.19), i.e., 

ipxx =Ac/?, (6.2a) 

</Jx =B(^, (6.2b) 

ipt = - Aipxxxj (6.2c) 
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and 

A = 11. (6.3) 

A is the auxiliary matrix that we introduced to deal with the complex operation in (6.2b) and it 
works in practice. As a result, with the help of A we solved the above CES and then categorized 
the solutions to the mKdV equation in terms of the canonical form of A (rather than canonical 
form of B). Solutions are categorized by solitons (together with their limit case) and breathers 
(together with their limit case). There are no rational solutions arising from (6.2) because no 
rational solutions correspond to zero eigenvalues of A while we need |B| 7^ to finish Wronskian 
verification. To derive rational solutions for the mKdV equation (1.3), we employed the Galilean 
transformed equation, i.e., the KdV-mKdV equation (1.4) which admits rational solutions in 
Wronskian form. Then the rational solutions to the mKdV equation can be recovered through 
the inverse transformation. Dynamics of some obtained solutions was analyzed and illustrated. 
Here, particularly, we would like to sum up a typical characteristic of limit solitons: the wave 
trajectories asymptotically follow logarithm curves (combined with linear functions). This point 
is based on several examples we have examined [43]. 

Obviously, through the Galilean transformation (4.2), all these obtained solutions of the 
mKdV equation (1.3) can easily be used for the KdV-mKdV equation (1.4) which often appears 
in physics contexts. In fact, in the paper we do not differ them from each other. In addition to 
the KdV-mKdV equation, our treatment to the complex operation in (6.2) can also be applied 
to the sine-Gordon equation. 

There are Miura transformations between the KdV equation and the mKdV equation (for 
both e = ±1). For the mKdV(e = —1) equation the Miura transformation provides a real map 
between solutions of the KdV equation and the mKdV equation(e = —1), c.f. [20]. However, 
when e = 1, i.e., for the mKdV equation (1.3), the Miura transformation has to be complex. In 
more detail, it maps the real mKdV equation (1.3) to a complex KdV equation. That means we 
have had a nice determinant expression for the complex KdV equation. Further investigation 
about this will be considered separately. 

6.2 List of solutions 

Let us list out the obtained solutions and their corresponding basic Wronskian vectors. Solutions 
to the mKdV equation (1.3) can be given by 

V = 2(arctan|)^ = "^^^'gf^'^f ^''^^ (6.4a) 

where 

/ = /(v?) = |iV^l|=Fi + iF2, Fi = Re[/], F2 = Im[/]. (6.4b) 

The available Wronskian vectors are the following. 

• For soliton solutions: 

V' = V^lv = ('^i''/'2,--- ^V'Af)^, (6.5a) 

with 

"Pj = 4^^' + ia]e-^^ , ij = kjx - Akp + if\ a+, aj,kj,cf G M. (6.5b) 
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• For limit solutions of solitons: 

^ = ip^l^\ki) = A-^Q+ + iA-Qo, ^^gGjv(M), (6.6a) 

with 

Q^ = (Qjo>2jir--,Qj^-if, Qt = l^9l^^^\ (6.6b) 

where ^i is defined in (6.5b). 

• For breather solutions: 

^ = V2N = {^11^^12, V521, V'22, • • • , ^N1,^N2V, (6.7a) 

with 

ipji = o,je + bj^^ ■, Vj2 = o-j^ — ^je~^^ , (6.7b) 

Cj = kjx - 4k]t + Cf\ aj, bj,cf G C. (6.7c) 

• For limit solutions of breathers: 

and the elements are given through 

9^" = {'Pi,2^'f2,2r-- ,'f']^^2)^ = AQ^ - BQ^ , (6.8c) 

where A ■Be Gn{C), 

Qt = (2jo> 2ji, • • • > QtN-if^ Qts = ^^fci^^^^ (6.8d) 

and ^1 is defined in (6.7c). 

We note that, thanks to the hnear property of the CES (2.18), one may also get mixed 
solutions by arbitrarily combining the above vectors to be a new Wronskian vector. For example, 
take 



^nA^n^- 



i) 



V2 

The related solution corresponds to the interaction between A'^i-soliton and a (A''2 — l)-order 
limit-soliton solutions. 

Finally, for the rational solution, it is given by 



2(Fi, xF2 - Fi F2,x) ^ _ ^ ^2. 
F| + Fl 



v{x, t) = VQ- ^ ' ' — ' ' % X = x- Qvlt, vo / G M, (6.10a) 



where still 

f = f{^) = \lf^l\=Fi+iF2, Fi = Re[/], F2 = Im[/], (6.10b) 
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and the Wronskian is composed by 

V' = (V'i,V'2,--- ,V'Af)^, (6.11a) 

with 

(i = 0,l,---,iV-l), (6.11b) 



V'i+l = TTTTTTTTT^V^l 



fci=0 



(2j)! aA;i2^' 
and 

ipi = \/2t>o + 2iA:i e''^ + ^2^0 - 2iA;i e"''S r?i = fciX - 4A;?t, /ci G M. (6.11c) 
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A Proof of Theorem 1 

Proof. The compatibility of (2.12a) and (2.12b), i.e., (/)^i = (/>te, yields (2.13b). Using (2.12a) 
one gets the complex conjugate form of / as 

/ = |^(t)||-l,]V^|. (A.l) 

Then the necessary derivatives of / and / are presented as the following, 

/,=|iV^,iV|, (A.2a) 

U^=\N^,N-l,N\ + \N^2,N + l\, (A.2b) 
U^^=\lf^A,N -2,N -l,N\ + \lf^2,N + 2\ 

+ 2|]V^,iV- l,iV + l|, (A.2c) 
f^ = - 4(|]V^, N -2,N -l,N\- \N~^?,, iV - 1, A^ + 1| 

+ \N^2,N + 2\) + iT{C{t))\N^ll (A.2d) 



and 



U=\B{t)\\-l,N-?,,N-ll (A.3a) 

/,, =\B{t)\{\ - 1, A^n, AT- 2,Ar- 1| + I - 1, aTTs, AT]), (A.3b) 

h^^ =\B{t)\{\ - I, N^, N + l\ + 2\- I, N^A,N- 2, N\ 

+ 1 -l,iV^,Af-3,Ar-2,Ar-l|), (A.3c) 

/i = - 4|5(t)|(| - 1,]V^, , AT - 3, AT - 2, AT - 1| 
-I -l,iV^,Af-2,Ar| + |-l,]V^,Ar + l|) 
+ \B{t)\t\ - l,N^2\+iT{C{t))\B{t)\\ - l,N^2\. (A.3d) 
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Using the condition (2.12a) the complex conjugate of (A. 2d) is 

/i = -4\B{t)\{\ - 1,N^5,N-3,N-2,N-1\ - \ -l,lf^A,N-2,N\ 

+ I - l,ff^3,N + l\)+tr{Cit))\Bit)\\ - l,iV^|, 

which should be same as (A. 4). This requires \B{t)\t = and tr(C(t)) G M(t), i.e., the condition 
(2.13a). 

Noting that (pxx = B{t)B{t)(/) and using Proposition 1 with Qj^s = d'^, we find 

tT{B{t)B{t))\ -1,N^3,N- 1| = -I -l,lf^5,N-3,N-2,N -l\ 

+ \-l,N^3,N + l\, (A.4) 

tT{B{t)B{t))\IV^,N\ = -\N~^A, N - 2, N - I, N\ + \N^,N + 2\, (A.5) 

tr{B{t)B{t))\ -l,lf^2\ = -\ -l,]\r^,iV-2,iV-l| + | -l,iV^,A^|, (A.6) 

tT{B{t)B{t))\lf^l\ = -\lf^3, N-l,N\ + \N^2, N + l\. (A.7) 

Then, substituting (A. 2) and (A. 3) into (2.10a) and making use of (A.7), we have 

It J J It ~T~ J XXX J 'J J XX J X ~r ojxjxx J Jxxx 

= 6\B{t)\{- I - l,lf^3,N + l\\lf^l\ - I - l,lf^A,N-2,N- l\\lf^2,N\ 

+ \N^A,N-2,N- 1,N\\ -l,lf^2\ - \N^3,N- l,iV + l|| -l,lf^2\ 

+ \N^2,N + 1\\ - l,lf^3,N- 1| - I - l,N^A,N-2,N\\lf^\), 

which is zero in the light of Proposition 2. Similarly, one can prove (2.10b). Thus the proof is 
completed. D 

B Eigen-polynomial of A = BB 

We prove the Proposition 3 through the following 2 Lemmas. 

Lemma B.l. For two arbitrary Nth-order complex matrices A and B, 

det(A/7v - AB) = det(A/jv - BA), (B.l) 

where I^ is the Nth- order unit matrix. 

Proof. Assuming that rank(A) = r, then there exist A^th-order non-singular matrices P and Q 
such that 

■4 = p(J»)e. (B.2) 



Thus 



det{\lN-AB) = det(A/jv-P( p" [J ] QB) 



.^U^T,. 1,1 Ir 



det{p-\XlN - P [ J Q ]QB)P) 



,-U^i Ir 



det(AIjv-P-^(P(^Q ^JQB)P) 
det(A/jv - ( Q*- 1)qBP). 
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In a similar way we have 

det(A/7v - BA) = det{XlN - QBP (^^ M ). (B.3) 

If we rewrite the matrix QBP into the foUowing block matrix form with same structure as 
Ir 


Qj^p = ( I'' f' ) , (B.4) 

V -D21 -D22 



then we have 



Ir 




)QBp^(f' „^'0.<3Bp(^' 2 )-(«:: o)' <"-^' 



which further means 

det(A/7v - AB) = A^-''det(A/r - Bu) = det(A/7v - BA). (B.6) 

We complete the proof. D 

Lemma B.2. Assuming that B is an arbitrary Nth-order complex matrix and B is its complex 
conjugate, then det{\I — BB) is a polynomial of A with real coefficients. 

Proof. Write 

/(A) = det(A/^ - BB) = onX^ + • • • + aiA + oq- (B.7) 

Then using Lemma B.l with A = B we have 



/(A) = det(A/^ - BB) = det(A/jv - BB) = det(A/7V - BB) = /(A), (B.8) 

which means all the coefficients {oj} are real. D 

C Discussions on the trivial solutions to the CES (2.18) 

Let the square matrices A,B follow the relation 

A = B]B. (C.l) 

We start from the 2x2 case. Noting that the product of all the eigenvalues of A is non- 
negative, in the following we first look at 

A = r ^^1 ^p y h^k2^0, ki,k2eR, (C.2) 

and suppose 
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with undetermined a,b,c,d£ C. However, in this case it can be found that the matrix relation 
(C.l) does not have any solutions unless kf = /cl- So next we turn to consider 

A=(^^' °^2 ), k€R. (C.4) 

In this case, the equation (C.l) admits a non-diagonal matrix solution B as 

e'\ c/0, c,d,eeR. (C.5) 



-d -=-^:f^ \ je 



c d 

We note that such a B does not lead to any nontrivial solutions to the mKdV equation. In fact, 
in the CES (2.18), the general solution to the equation set (2.18a) and (2.18c) is 

with arbitrary matrix H € C2x2 and 

i = kx + Akh + ^^^\ A:,e^°)GM. (C.6b) 

However, no matter what condition the matrix H should satisfy under the equation (2.18b), the 
Wronskian 

is always a constant, which leads to a trivial solution to the mKdV equation. In the case of the 
N X N matrix (N is even) 

^ = diag(-A;^-A:^••• ,-fc^). A; G M, 

similar to (C.6), the general solution to (2.18a) and (2.18c) is 



H x{e'^,e-'^,0,0,--- ,0f, HeCN 



xN, 



which leads to a zero Wronskian f{<f)- 

With these discussions we can conclude that for the CES (2.18) with an A^ x A^ matrix 

A = diag(a,a, • • • ,a), a € M, (C.7) 

the possible solution (p to (2.18) composes a trivial Wronskian f{(p)- 

D Proof for Theorem 2 

Proof. The parameter vq will lead to complicated expressions for / and its derivatives. For 
simplification let us introduce the notation | • \j where the subscript j indicates the absence of 
the <p^^' column[22], for example. 
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Derivatives of / have already given in (A. 2). For /, using the CES (4.6) we can reach 

N 

f=\B-Ht)\Y,{-voyi''-^\N\„ 

j=0 

N-l N-1 

j=0 j=0 

N-l 

fxx =\B'\t)\{{-vo)''\N^2, AT + 1| + ^ ^_v^yiN-j\ffri^ ^ + 2\j 

j=0 

N-2 N-2 

+ Y. {-voyi''-'\N^2,N,N + 1|, + Y {-voy^^i''-'-^\N\j 

j=0 j=0 

N~2 

+ 2Yi-voy^'i''~'~'\N^l,N + l\j), 

j=0 

N-3 

fxxx=\B-\t)\{3Yi-^oy^^^''-'-^\N^,N + l\, + i-vof\N^,N + 2\ 

j=0 

N-2 

+ 3 J^ {-voy^H^~^~^\N^l, N + 2\j + 2(-t;o)^|iV^, N-1,N + 1\ 

j=0 

N-l N-2 

+ Yl i-voyi^'^\N^l,N + 3\j + 2Y {-voyi^-^\N^2,N,N + 2\j 

j=0 j=0 

N-3 

+ 3 ^ {-voy+^i^-^-^\N^, N,N + l\j + 2{-vo)^-^i\N^3, A^, iV + 1| 
i=o 

N-3 N-3 

+ Y {-voyi^-'\N^,N-l,N,N + l\, + Y {-voy^^i^-^-^\N\j), 

j=0 j=0 

N-l N-3 

ft = -4\B-\t)\{Y {-voyt^'^'iN^hN+si, + Y i-voy^^^-'-^\N\j 

j=0 j=0 

N-3 

+ Y {-'^oyi^~^\N^3, N-l,N,N + l\j + {-vo)^\N^2, N + 2\ 

j=0 

N-2 

- Y {-voyi^'^\N^, N,N + 2\j - i-vo)^ '^i\N^3, iV, iV + 1| 

N 

- {-vof\N^3, N-l,N + l\)+ tr(C(t))|S-i(t)| Y {-voyi''~'\N\j. 

3=0 
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Using the condition (4.6) the complex conjugate of (A. 2d) should be same as the above ft- This 
requires tr(C(t)) G M(t), i.e., the condition (4.7b). 

Besides, noting that (j)xx = {B{t)B{t) — VqIn)(I> &iid using Proposition 1 with Qj^s = d'j^, we 
can have the following identities: 

tT{B{t)B{t) - v'^In)\N\j =|iV^, N + 2\j - \N^2, N,N + l\j 

-\N\j-2, (i = 0,l,---,iV-2), 

tr{B{t)B{t) - v^In)\N^I, N + 1\j =\N^l, N + 3\j - \N^3, N-l,N,N + l\j 

-|iV^l,iV + l|,-_2, (i = 0,l,--- ,iV-3), 
tT{B{t)B{t) - v'^In)\N^I\ = - \N^3, iV - 1, iV| + \lf^2, N + l\, 
tj:{B{t)B{t) - vilN)\N^, N\=- \N^, N-2,N-1,N\ + \N^2, N + 2\. 

With these results and (A. 2) in hand, for (4.4a) we have 

ftf - fft + fxxxf - 'ifxxfx + 3/x/xx - ffxxx 

=6{-vo)^~^ii\lf^3, N, N + l\\lf^l\ - \lf^2,N\\N^3,N- l,iV + 1| 

+ |]V^,iV + l||]V^,iV- l,N\) 
=0. 

Similarly, one can prove (4.4b). D 
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